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Abstract 

We study CP violation in a multi-Higgs doublet model based on a 5*3 x Z3 
horizontal symmetry where CKM phase is not the principal source of CP 

m 

violation. We consider two mechanisms for CP violation in this model: a) 
CL|| CP violation due to complex Yukawa couplings; and b) CP violation due to 

Q ' scalar-pseudoscalar Higgs boson mixings. Both mechanisms can explain the 

observed CP violation in the neutral Kaon system, e'/ e due to neutral Higgs 

•i— I , 

^ ' boson exchange is small in both mechanisms, but charged Higgs boson con- 

^' 3 . 

tributions can be as large as 10 for a), and 10 for b). CP violation in 

the neutral B system is, however, quite different from the Minimal Standard 
Model. The neutron Electric Dipole Moment can be as large as the present ex- 
perimental bound, and can be used to constrain charged Higgs boson masses. 
The electron EDM is one order of magnitude below the experimental bound 
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in case b) and smaller in case a). 
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I. INTRODUCTION 



The origin of CP violation, fermion masses and mixings are some of the outstanding 
problems of particle physics today. The Minimal Standard Model (MSM) provides no ex- 
planations for the observed fermion masses and mixings. It is believed that one needs to go 
beyond the MSM to solve these problems. There is no satisfactary solution so far. The best 
one can do at the present is to reduce the number of free parameters in the theory regarding 
fermion masses and mixings. The understanding of CP violation is also very poor. So far 
CP violation has only been observed in the neutral kaon system. In the MSM CP violation 
is from the phase in the CKM matrix Vkm- The model is consistent with observations. 
It is, however, important to study different processes and other models of CP violation to 
better understand the origin of CP violation. In this paper we will study some details of 
CP violation in a model proposed by Ma 0. The model is based on a horizontal 5*3 x 
symmetry. In this model there are some interesting relations between the fermion masses 
and mixings. The free parameters related to the fermion masses and mixings are reduced. 
This model also has some interesting consequences for CP violation 0,^. Many of the 
considerations in this model will have applicability to more general multi-Higgs models. 

There are four Higgs doublets in the x Z3 model. Their interactions at the tree level 
mediate neutral flavor changing current. Like any multi-Higgs doublet models with flavor 
changing neutral current at the tree level, there are different mechanisms for CP violation. 
CP violation can arise in three places in this type of models: 1) Non-trivial phase in the 
Vkm matrix; 2)Non-trivial phases in the flavor changing Yukawa couplings; and 3) Mixings 
of scalar and pseudoscalar Higgs bosons. In cases 2) and 3), CP violation can occurs at 
the tree level by exchanging neutral Higgs bosons These models have much richer 
phenomenology for CP violation than the MSM. In the MSM, CP violation in the neutral 
Kaon system can be explained by exchanging W-bosons at the one loop level (the "box 
diagram") 0. There are specific predictions for e'/e resulting from the direct CP violation 
in Kls ~^ 27r [0; and predictions for CP violation in the neutral B system The Electric 
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Dipole Moment of neutron is generated at the three loop level, and is less than 10~^^ecm 
0. The electron EDM is even smaller |T0[. In the 5*3 x model, the CP violation coming 
from the phase in the CKM matrix is inadequate to account for the e parameter in the K 
system. One of the other two has to be invoked. We consider the consequences of either 
mechanism for: (i) e'/e; (ii) CP violation in the neutral B system; and (iii) the neutron and 
electron EDMs. The results are dramatically different from the MSM. 

II. YUKAWA COUPLINGS IN THE S3 x Z3 MODEL 

In the 5*3 x Z3 model, there are four Higgs doublets, 0i,2,3,4- The quarks and Higgs bosons 
transform under the 6*3 x symmetry as [0] 

QsL ,tR,bR,4>i: (1, 1) , 
(giL,g2L) , (03,04) : (2,u;) , 
{cr,Ur) ,SR,dR) : {2,uj^) , 

02:(1,^'), (1) 

where 7^ 1, cj^ = 1 is the Z^ element. The Yukawa couplings consistent with the S3 x Z3 
symmetry are given by 

Ly = -fl{qiL<p3UR + g2L04Cij) - /2g3L01^i? 

- /3(giL02SR + q2L<p2dR) - hiqiL^bR + g2L04&i?,) 

- fbiqsLhdR + qzL(t)iSR) - f6q3L(pibR + H.C. (2) 

where (pi = (0°*, —(p^)'^. Without loss of generality we work in a basis where all Vaccum 
Expectation Values (VEV) are real. When the neutral components develop VEVs, < (pi >= 
Vi, we obtain the quark mass matrices 

ffivs \ 



fm 
V f2vj 
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The quark phases can be chosen such that 
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with a, b, c, d real and ^ = l^le*"^ complex. M'^ can be diagonalized by a bi-unitary 
transformation 



(5) 



Here M'^ is the diagonalized down quark mass matrix. Vl and Vr are unitary matrices. 
Because the up quark mass matrix is already diagonlized, Vl is the CKM matrix Vk hi- 
lt is convenient to work in a basis of the Higgs bosons in which the Goldstone bosons 
are removed. To this end we define the following 
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H2 



(6) 



1)124 ™i24 

where = vf + v^, ff24 = "^i + "^i + "^45 f ^ = + f| + f| + f|. The transformation 
is the same for both the neutral and charged Higgs bosons. For the neutral Higgs bosons, 
G = + iGzi where Gz is the Goldstone boson 'eaten' by Z, and is a physical field 
whose couplings are the same as the Higgs boson in the MSM. For the charged Higgs bosons 
G is the Goldstone boson 'eaten' by W. In this basis, we have 

Dpl,0 



- ULVKMYfDRhf + DLVlj.jYl'URh- + H.G. 



(7) 



where hi are the component fields of with — [hf ,hy \/2). Ul,r 
Dl,r = {d, s, b)J^j^. The Yukawa couphngs are given by 



VuVl 



Y2 = Diag{0, , ) , 
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In general /i^'"*" are not the mass eigenstates. They will mix with each other. In particular 
if CP is violated in the Higgs potential, i?e/i° and /m/i° will mix. Also the charged Higgs 
boson mixing matrix will be complex. We can parametrize the mixings as 

/ /i° \ / aoo I3qj \ f rio\ 

ftfeO ttfcj P'kj 



Rehl 



\Imhl) 



\Ij ) 



(9) 



where and r]i are the mass eigenstates, the matrix {a(5) is a 7 x 7 othogonal matrix, 

and (7) is a 3 X 3 unitary matrix. 
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The specific numbers for aij, a'^j, j3ij and depend on tlie details of the Higgs potential 
Unfortunately they are not determined. To simplify the problem, we will discuss two cases: 
a) CP violation only come from complex Yukawa couplings; and b) CP violation only come 
from the mixings of real and imaginary /i?. Case a) can be realised by constraining certain 
soft symmetry breaking terms in the potential [Q. We further assume, for simplicity, that 
Reh^^ are the mass eigenstate Ri and consider their effects. The same analysis can be easily 
carried out for /m/i° in the same way. The source for CP violation is the non-zero value 
for cr which is a free parameter. We will present our results for a = 80", which is close to 
the maximum of the allowed phase. Case b) can be realised by requiring spontaneous CP 
violation. The value of a will be zero and CP violation arises due to scalar-pesudoscalar 
Higgs boson mixing. For illustration, we consider the effects of a neutral mixed state 

R = cosORehl + sinOImhl , (10) 

and for the charged Higgs boson we consider mixing 

r]+ = -f22ht + i2'iht (11) 

where are complex numbers, and I722P + I723P = 1- We assume that R and r]'^ are mass 
eigenstates, and all other Higgs bosons are much heavier and their effects can be neglected. 

The details of the results depend on the specific vules of a in case a) and how Higgs 
bosons are mixed in case b). However, the special cases considered here will serve as a guide 
for a more complete analysis if the details of the mixings are known. The general features 
will be the same. We will comment on other cases later. 

The parameters a, 6, c, and d are constrained from the down quark masses and the 
CKM mixings. We take as input parameters a = O.OAGeV, b = 0.25GeV, c = 2.6QGeV , 
d = 4GeV. The mass eigenvalues for the down quarks are quite insensitive to the phase a. 
For both cases, we have = A.SGeV, rris = lAdMeV and = 9.5MeV. These values are 



well within the allowed regions WU. The CKM matrix for case a) is 
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/ 0.975 -0.222 0.00476 \ 

Vkm= 0.221 +i0.0033 0.974 + i0.014 0.043 - iO.0015 , (12) 

V -0.014 + il. 2 X 10-5 -0.041 -z6.8 X 10-^ 0.998 -^0.034 / 
and for case b) 

/ 0.975 0.22 0.0048 \ 
Vkm = -0.219 0.975 -0.0436 . (13) 
V -0.014 0.0415 0.999 j 
The values for the VEV's are not fixed, we only know vs/v^ = rriu/mc. We will use the 
values in Ref [||: f i = f 2 = MGeV, = 0.9 and V4 = 238GeV for illustration. We shall 
comment on effects of changing these values later. 



III. CONSTRAINTS ON THE HIGGS BOSON MASSES FROM THE NEUTRAL 

K AND B MESON SYSTEMS 



The 5*3 X Z3 model has very restrictive allowed values for the non-trivial CP violating 
phase in the CKM matrix. The CP violating measure J [T^ is less than 2.5 x 10"^. If CP 
violation is purely from the CKM matrix, it is not possible to explain the observed CP viola- 
tion in the neutral K system which requires J > 10^^. Therefore in this model CP violation 
due to Higgs boson exchange has to be considered. Because the neutral Higgs bosons medi- 
ate flavor changing neutral current at the tree level, there are stringent constraints on their 
masses and their interactions arising from experimental data. Some of the best constraints 
are from the mass differences in the neutral K and B systems. We must make sure that the 
observed CP violation in the neutral K system e = 2.3 x lO^^e*'^/^ is explained. 

The CP violating parameter e is given by 



'12 ^t7r/4 



(14) 



V2AmK 

where is the matrix element which mixes with K^, and Aitlk is the mass difference 
between and rriKg- The AS* = 2 Hamiltonian, responsible for M^, generated by 

exchanging neutral Higgs bosons Ri is given by 
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Heff = (d[{a,. + tai,)Y,%)^^ + - . (15) 



We obtain 



= < K°\Heff\K^ > 

flmx f 5 m\ |., , , - ^ -2 I / -I \i/d* ^2l 



+ (a,, + ^a',,)Y,%{ay, - ia',,,)Y,%,{- + ., )) . (16) 

Here we have used the vaccum saturation and factorization approximation results for the 
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— — — 1 ?T7^ 

< K'\d{l±^,)sd{lTl,)s\K' > =fl:mK{-+ ^ 



^6 {ms + rridY 

< K'\d{l ± 75)3^(1 ± j,)s\K' > = -IflniK , ,2 ■ (17) 

6 (ms + nidY 

The contribution to the mass difference /\mK is given by 2ReMi2. Similar formula holds 
for the neutral B system. 

To constrain the Higgs boson masses, we require that the neutral Higgs boson contribu- 
tions to the mass differences in the neutral K and B systems to be less than the experimental 
values: Amx/niK = 7 x 10~^^, and AmB/mB = 8 x 10~^^. We find that for case a) the 
tightest constraints on the masses of Reh^ 2 are from the mass difference AMb of the neu- 
tral B mesons which gives Mh^ > 2.9TeV and Mh^ > S.lTeV. With these masses, Reh^2 
can not produce large enough e. Similar consideration yields Mh3 > 3.5TeV, and we find 
the experimental value of e can now be produced if the mass is about 5.6TeV. The mass 
difference AMk of the neutral K mesons gives weaker bounds in all cases. For case b), the 
experimental value of AMb constrains Mji > 3TeV. From the experimental value of e, we 
obtain sinOcosO / M^^ = 1.1 x 10"^ which implies Mr < 7TeV. 

Similar analyses for case a) and case b) have been carried out previously in Ref. 0] and 
Ref . 0] respectively. Our analysis for case a) is similar to that in Ref. @] . We used different 
values for the paramters a, b, c, and d which obtain smaller mass for rrih (A.SGeV) compared 
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with the value bAGeV used in Ref. 0]. The analysis for case b) is different than that in 
Ref. 0. In Ref. [Q, the averaged effects of Higgs boson exchange were considered. 

From the above we see that the Higgs boson masses are constrained to be in the multi 
TeV region. One would wonder if such heavy Higgs bosons may violate the unitarity bound. 
However this is not the case. It should be noted that the unitarity bound only apply to 
Higgs boson mass. Its mass is not constrained in the cases we are considering, can be light. 
In the above discussions we have neglected mixings between with h^. If such mixings are 
large, the unitarity bound can constrain the masses. However there is enough freedom 
in our model to make the mixings with Hq sufficiently small such that the unitarity bounds 
are always satisfied This argument applies to many models. A well known example is 
the minimal supersymmetric standard model. In this model, when the soft supersymmetry 
breaking parameter goes to infinite, the other heavy Higgs bosons decouple from the theory 
and cause no problem to the unitarity bound n . 



IV. PREDICTIONS FOR e'/e. 

In this section we study the direct CP violation in Kl s ^ 27r decays. CP violation in 
these processes is characterized by the value of e'/e. e'/e is defined as 

e^ _ - ImA2/ReAo 

e ~ ^ + ItuMu/AMk ' ^ ^ 

where u = ReA2/ReAo = 1/20, ^ = ImAo/ReAo. Here Aq and A2 are the AI = 1/2, 3/2 
decay amplitudes for K^ s ~^ Svr. 

In our model, the neutral Higgs boson can induce ImAi at the tree level. These am- 
plitudes are constrained to be very small due to large Higgs boson masses. They are small 
compared with the CP conserving amplitudes ReAi generated by W-boson exchange at the 
tree level. The neutral Higgs boson contributions to e'/e are very small. However there may 
be large contributions from the charged Higgs bosons. The charged Higgs boson contribu- 
tions to ImAi can be generated at the one loop level. The dominant one is from the charged 
Higgs boson mediated gluon penguin. The relevant AS* = 1 effective Larangian is given by 
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Las=i = tda^'ifi^^ + /2^^)A^G;:, , (19) 



where is the gluon field strength, and 

.2 



/i+ h 



3 

/2 = - ^^^)^"^{(^^MV;S.)r2(n"VK./7.,)n}C/ , (20) 



where Q = (asimh) / as{iJ^)Y'^^'^^ ~ 0.17 is the QCD correction factor, and 1 is summed over 
u, c and t. We will use as{fi) ~ An/G for fi = IGeV. The above effective Lagrangian will 
generate a non-zero value for ImAo Las=i also generates a non-zero value bld for CP 
violation in and mixing due to long distance interactions through and vr, rj, rj' 



mixings [0. One obtains [|T6|,|T7| 



2mKlmM^^^Lj, ^ 0.8 x IQ-\J\ + /2)(Ger3) , (21) 

where D is a supression factor of order 0(m|^, m^)/A^. ^/ein is of order -0.014 to -0.1. 

We find that in both a) and b) cases, the donimant contributions are from the top quark 
in the loop arising from mixing in the charged Higgs boson couplings. For case a), we have 



, GeV"^ mt , m? 



iMht) « -T^T^ta^ . (22) 

And for case b), we have 



Ho « -7.35 X mM-.^ra^j^ln^ • (23) 



The contributions to e can be significant in both cases depending on the Higgs boson masses 
and the CP violating parameter 1771(722723) • If the masses of the charged Higgs bosons are 
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of order a few hundred GeV, Ild can be as large as the experimental value and e'/e can be 
of order 10~^. However, there are constraints on the charged Higgs boson masses from the 
experimental upper bound on the neutron EDM. We will study these constraints in Sec. VI. 
When these constraints are taken into account, eiD is generally constrained to be less than 
3 X 10^^ for case a), e'/e is then constrained to be less than 3 x 10~^. However, for case b), 
eLD can still be as large as 10"'^ and e'/e can be 10^^. 

V. CP VIOLATION IN THE NEUTRAL B SYSTEM. 

There are many processes which can test CP violation in the neutral B system. Some 
particularly interesting ones are p] 

B, ^ J/^Ks , 

Bd ^ vr+TT- , (24) 
Bs ^ pKs . 

The differences of time variation of decay rates for the above processes and their CP tran- 
formed states are given by 

TiB%t) ^ fcp) - T{B%t) fcp) 



O'fCP 



TiB%t)^fcp) + nB%t)^fcp) 



(1- 




\'^)cos{AMBt) 




2ImXsin{AMBt) 


1 + 


\M 


2 



(25) 



where fcp indicates the final states. A is defined as 



where {q/p)B = \J / My^^ A and A are the decay amplitudes. If the final state contains 
Ks, S = {q/p)K which has a phase of order 10"'^. For other cases S is equal to one. 

Non-zero asymmetry a/cp signals CP violation. If |A| is not equal to one, it indicates 
that CP is violated in the decay amplitudes. In the MSM |A| is equal to one to a very good 
approximation for the above three processes. The asymmetries are proportional to ImX. 
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In the MSM, the processes in Eq.(24) measure three angles related to CP violation in the 
CKM matrix, 



ImX{Bd J/i^Ks) = -sin2(3 , 
ImX{B(i Tr"*"?!") = sin2a , 
Im\{Bs — > pKs) — —sin2^ , 



(27) 



where 



V, 



KM,tdVKM,tb 
VKM,udyKM,ub/ 
VKM,cdVKM,cb\ 



a — arg 

(3 = arg 
7 = arg 



KM,tdyKM,th J 
VKM,udVKM,ub 



(28) 



, VKM,cdVKM,cb) J 

The sum of these three angles is equal to 180^. 

In the 5*3 x model, the situation is very different. In this model CP violation is 
mainly due to neutral Higgs boson exchange. The CP violating decay amphtudes A and A 
are small because the decay amplitudes are dominated by the CP conserving tree level W 



interactions. However the phase of ^Mj^*/M]| in the B — B mixing due to neutral Higgs 
boson exchange can be large. In case a), there is CP violation arising from the phase in 
Yukawa coupling of Higgs bosons, as well as CKM matrix, but the former is much larger. 
The three meaurements in Eq.(24) do not measure the angles a, (3 and 7 defined in Eq.(28) 
anymore. The first two processes will mostly measure the phases in M-^^ . We have 



Im\{Bd TT+TT") fa ImX{Bd J/^Ks) < 0.42 , from Rehl, 
Im\{Bd 7r+7r~) f« ImX{Bd J/^Kg) < 0.19 , from Rehl , 
ImX{Bd TT+TT") fa ImX{Bd J/i^Kg) fa 0.19 , form Rehl ■ 



(29) 



For case b), the CP violation is purely from scalar-pesudoscalar mixing in the Higgs 
sector and we find 



ImX{Bd Tr+TT") fa ImX{Bd J/i^Kg) fa -0.25 



(30) 
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ImX for Bs pKs is different for a) and b). For case a), the neutral Higgs boson 
contributions to tlie asymmetry are small. However ImX{Bs pKs) due to CP violation 
in the KM-matrix can be about 0.1. For case b), Im\{Bs — > pKs) from neutral Higgs boson 
exchange is only about 0.02. 

If interpreted as in Eq.(27), we find for case a), sin2a = —sin2P, sin'y = 0.05, and 
a + /? + 7 7^ ISO''. For case b), we have, sin2a = — sm2/?, siri'y = 0.01. We again find, 
a + /5 + 7 7^ ISO''. These predictions are different than those of the MSM. 



VI. THE NEUTRON ELECTRIC DIPOLE MOMENT. 

The prediction for the neutron EDM in the S3 x model is very different from the MSM. 
In both a) and b) cases, the neutron EDM can be generated at the one loop level, and the 
two loop contributions can be significant. The one loop contributions in this model are 
also different from multi-Higgs models with neutral flavor conservation for tree level neutral 



Higgs boson exchange [0. In flavor conserving models the fermion EDMs are proportional 
to the third powers of the external fermion masses. For u- and d- quarks and electron, the 
EDMs at the one loop level are very small (< lO^^^ecm). In the models considered by us, 
the couplings are quite different from the flavor conserving models. Furthermore the off 
diagonal couplings may contribute significantly. It is possible to have a large neutron EDM 
at the one loop level. The d quark EDM due to neutral Higgs boson exchange at the one 
loop level is given by 



Qde ^ , . mi 3 (m^ 

--—rIm{ai diai id)^-{- - In \ — 
16n^ mi, 2 \ 



dd = ^7-^Im{ai^dia^d)^^Q - In ( '-^ ))Cd , (31) 



Ri 

where ai^gi = {aki + ia'i^i)Y^gi, (d = («s('^h)/«s(Ai))^^''^'^ ~ 0.12, and 1 is summed over d, s, 
and b. The quark EDMs are related to the neutron EDM by quark model, 

Dn = Md/3 - du/3 . (32) 

For case a) we have 
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Dn{d, Rehl) < 0.7 x lO^^^ecm , (33) 
Dn{d, Rehl) ^ 2 x lO^^^ecm . 

For case b), we have Dn{d) ^ 2 x 10~^^ecm. These values are at least three orders of 
magnitude smaller than the experimental upper bound of 1.2 x 10^^^ ecm ||2T|. The u quark 
EDM is zero at the one loop level. 

In multi-Higgs models, there may be large contributions to the neutron EDM at the two 
loop level from the Weinberg operator |]18| Dn{W) and from the color dipole moment of 



gluon due to Bar-Zee type of diagrams [T^j2^ Dn{BZ). In our model, we have 



1 2 2 

Ttl Tfi 



647r2 "m^o ml, ' 



where Cvk ~ 6 x 10^^, and Cfez ~ 10^^ are the QCD correction factors, Cu = 2 and Cd = 4, 
and A ^ IGeV is the chiral symmetry breaking scale. The parameters ImZ are given by 

ImZlt = -^ReiY^ssiaki - ia',,.^)Im{Y^,^^^{ak'i - ia'fc/J) , 
m^ 

I'rnZl^ = Jm(y;^33(afci - m^Jl^^ ^("fc'i " ^a'fc'J) 

rriumt 

I ~ 

ImZla = Im{Y,^^^{ak^ - ia'k^)yk'^Mk'^ + ia'k'i)) ■ (35) 

mdmt 

For case a), because there is no CP violation in the up quark sector only down quark loops 
contribute, Dn{W) from the Weinberg operator at the two loop level is small. There are 
non-zero Dn{BZ) from d-quark due to Bar-Zee mechanism. We find that the contributions 
from Rehi2 is also small (< 4 x 10"^^ ecm). Reh^ contribution is even smaller (< 10~^^). 

For case b), the two loop contributions to the EDM are significantly larger because in 
this case there is CP violation in the top quark interaction. We have 

DrXBZ,u) ^ (2 ~ 8) X lO^^^ecm , 

Dr,{BZ, d) ^ (2 ~ 8) X 10-27ecm , (36) 
15 



for rrit between 100 GeV to 200 GeV. The neutron EDM can be as large as the experimental 
upper bound. The contribution from the Weinberg operator is small, Dn{W) < 10~^°ecm. 

The charged Higgs bosons can also contribute to the neutron EDM. At the one loop 
level, the u and d quark EDM are given by 

1 771/1 TTV^ ~ -t- " 

'^'^ ^ i£^^^£^^^[7.-^7^^(^^^^')'i(^^Mn")u] . (37) 

Z^tTT 777', -L 777-, -l 
h. h. 

For du, 1 is summed over d, s, and b; and for da, 1 is summed over u, c, and t. At the two 
loop level, there is a large contribution from the Weinberg operator, 

Dn{W) ^ eCl^M^^ImZ'^^ln^ , (38) 

where ^(^ = 3 x 10~^ is the QCD correction factor, and 

ImZ'^ = -^Im[^^,^l{VKMYf),,{VKMY:),,] . (39) 

We find that in case a) the dominant contributions are from the two loop Weinberg 
operator. We have 

D^{W) « 1.6 X ecm , 



i9GeV\ ml ml 



Dr,{W) ^ 1.4 X 10-i«^^/n^— — ^ ecm , 

mi+ mi+ (ISOGeV/)^ 

D^{W) ^ 1.2 X ecm . (40) 

m^+ {15QGeV)^ 

Requiring the contributions to be less than the experimental value, we find the masses of 
hi 2 have to be larger than 2.5TeV. There is no constraint on mass. Combining this in- 
formation with those from Eqs.(22) and (23), we find the charged Higgs boson contributions 
to €ld is less than 3 x 10~^, and e'/e is less than 3 x 10"^. 

For case b), we find the dominant contribution is from the one loop d quark EDM. We 
have 

16 



DM ^ 5.4 X 10-^^ Im{^,,^;,)^ln^-^ ecm . (41) 



Requiring Dn{d) to be less than the experimental value, Eld is constrained to be less than 
10~^, and e'/e can still be of order 10~^. Assuming maximum mixing, the mass of r)~^ is 
constrained to be larger than 5 TeV. 

We also checked cases for different values of a and different mixings between Rehi and 
Imhj. We find that for case a) if 95° > \a\ > 85°, it is not possible to produce the 
experimental value for e because the constraints on the neutral Higgs boson masses from the 
mass difference of the neutral B mesons is too strong. The model can explain the observed 
CP violation in the neutral K system even for \a\ = 10°, and |180° — cr| = 10°. The neutral 
Higgs boson masses are typically constrained to be larger than a few TeV. For angles close to 
zero and 180°, CP violation for B system become small. ImX{Bd J/ipKs) ~ Im\{Bd 
7r"'"7r~) are between 0.05 ~ 0.45. Im{Bs —>■ pKg) is less than 0.1. The neutron EDM due to 
the neutral Higgs bosons is typically less than 10~^^ ecm. The contributions to the neutron 
EDM from the charged Higgs bosons can be close to the experimental upper bound, e'/e 
are typically less than 10~^. For case b), /mA(i?d — > tt+tt") and Im\{Bd J/ipKs) vary 
between 0.02 to 0.3. ImX{Bs — > pKs) is less than 0.1. The neutron EDM due to neutral 
Higgs bosons is larger than 10~^^edm and can be as large as the experimental upper bound. 
The charged Higgs bosons contribution to the neutron EDM can be close to the experimental 
upper bound, e'/e can be of order 10~^. 

The predictions also depend on the choices of the VEVs. The general features are, 
however, the same. For example for vi = 210 GeV, V2 = 3; GeV, = 0.5 GeV and 
V4 — 130 GeV, we find: l)The neutral Higgs bosons are constrained to be in the multi TeV 
region; 2) e'/e is small in case a) and can be 10~^ in case b); 3) the predictions for CP 
violation B system maintain the same features; and 4) the predictions for the neutral EDM 
from the neutral Higgs bosons are in the same regions as discussed earlier. 
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VII. THE ELECTRON ELECTRIC DIPOLE MOMENT 



The 5*3 X Z3 model may also have interesting CP violating signatures in the lepton 
sector. We assume the same 5*3 x Z3 assignments for the left handed and the charged right 
handed leptons as their quark partnenrs 0. The mass matrix and Yukawa couplings for 
the charged leptons are similar to the down quarks. One simply changes the parameters 
(a, b, c, d , and ^) for quarks to (ai ,bi ,Ci ,di , and = \^\e'^°'') for leptons. We use 0: 
ai = O.lOGGeV , 6, = , q = 1.781Gel^ , di = 8.6 x IQ-^GeV. For this set of parameters, we 
have rrie = 0.511MeV , = lOGMeV and m,- = 1784MeV which are in good agreement 
with experimental data. We again consider two ) CP violation purely due to complex 

Yukawa coupling with the phase a' = 80°; and case b) CP violation due to Reh2 and Imh\ 
mixing. If right handed neutrinos exist, the charged Higgs bosons will also contribute to 
electron EDM. However due to very small neutrino masses, the contribution to the electron 
EDM is neglegiblly small. We will not consider them here. 

For case a) we find that the one loop contributions are small (< lO^^ecm). However the 



two loop contribution due to Bar-Zee mechanism |19,22| can be as large as 



de{Rehl) < 1.5 x 10"^^ (42) 

for rrit < 2QQGeV . Reh^ contribution is much smaller. For case b), we find that the one 
loop and two loop contributions are small (< lO^^^ecm). 

For case a) the predictions depend on the phase a'. However as long as the phase 
is not too close to 0° and 180° (Isma'l > 0.17), the electron EDM can be a few times of 
10"^'' ecm. It is below the experimental upper bound p3[. For case b), different mixings may 
have different values for electron EDM and can be larger than the special case considered 
here, but is less than 10~^^ ecm. For example, the case for i?e° and Imh2, d^. is about 
2 X 10^^^ ecm. Varying VEVs can change the predictions for the electron EDM. Using the 
set of VEVs discussed at the end of the last section, we find the value for electron EDM is 
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smaller by one order of magnitude for case a) , and for case b) we again obtain small electron 
EDM. 



VIII. CONCLUSIONS 

We have studied in detail some effects due to two different CP violating mechanisms 
in the S3 x Z3 model. Both mechanisms discussed in this paper can explain the observed 
CP violation in the neutral K system. CP violation in the neutral K system and the mass 
difference in the neutral B system constrain the neutral Higgs boson masses to be in the 
multi TeV region. The predictions for other CP violations observables are very different 
from the MSM. 

i) e'/e : In the MSM, depending on the top quark mass, e'/e can be as large as 10~^. 
In the 5*3 x Z3 model e'/e due to neutral Higgs bosons is small. The charged Higgs boson 
cantribution for e'/e in case a) is also small. But for case b), e'/e can be as large as 10~^. 
The measurement of e'/ e may distiguish the MSM and case b) from case a) . 

ii) CP violation in B system: In section V, we discussed three asymmetries in the neutral 
B system. In the MSM model these three asymmetries measure three angles related to the 
CKM matrix. The sum of these angles is equal to 180°. However in the S3 x Z3 model, CP 
violation in the neutral B system comes from different sources. The predictions are very 
different from the MSM. If we interpret the measurements of the three parameters in Eq. 
(27) in terms of the three angles, we find their sums are not equal to 180° in both a) and 
b) cases and more over sin2a — —sin2p. These experiments should distinguish the MSM 
from the S3 x Z3 model and the B factory will provide us with very useful information. 

iii) The Electric Dipole Moment: The predictions for the neutron and electron EDMs 
are several orders of magnitude larger than those from the MSM. There are also dramatic 
differences between the two cases considered. The neutron EDM can be as large as the 
experimental upper bound for case b). It is smaller in case a). The electron EDM is below 
the experimental bound, but future experiments will reach the sensitivity necessary to test 
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case a). The electron EDM is smaller in case b). It is also interesting to note that some 
of the charged Higgs boson masses are also constrained to be in the TeV region. This is 
different from flavor conserving multi-Higgs models, where limits are much weaker. 
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